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Effects of interchain couplings and anisotropy on a Haldane chain have been investigated by single 
crystal inelastic neutron scattering and density functional theory (DFT) calculations on the model 
compound SrNi2V208. Significant effects on low energy excitation spectra are found where the 
Haldane gap (Aq « 0.41J; where J is the intrachain exchange interaction) is replaced by three energy 
minima at different antiferromagnetic zone centers due to the complex interchain couplings. Further, 
the triplet states are split into two branches by single-ion anisotropy. Quantitative information 
on the intrachain and interchain interactions as well as on the single-ion anisotropy are obtained 
from the analyses of the neutron scattering spectra by the random phase approximation (RPA) 
method. The presence of multiple competing interchain interactions is found from the analysis of 
the experimental spectra and is also confirmed by the DFT calculations. The interchain interactions 
are two orders of magnitude weaker than the nearest-neighbour intrachain interaction J = 8.7 meV. 

The DFT calculations reveal that the dominant intrachain nearest-neighbor interaction occurs via 
nontrivial extended superexchange pathways Ni-O-V-O-Ni involving the empty d orbital of V ions. 

The present single crystal study also allows us to correctly position SrNi 2 V 2 08 in the theoretical 
D-J± phase diagram [T. Sakai and M. Takahashi, Phys. Rev. B 42, 4537 (1990)] showing where it 
lies within the spin-liquid phase. 

PACS numbers: 75.50.Ee, 75.40.Gb, 75.50.Mm, 71.20.-b 


I. INTRODUCTION 

Spin-1 Heisenberg antiferromagnetic (AFM) chains or 
Haldane chains are of current interest due to their novel 
magnetic properties [1, 2]. In such systems, the mag¬ 
netic ions interact with their nearest neighbors in only 
one direction (ID) and the long-range magnetic order is 
suppressed even at T = 0 K by strong quantum spin 
fluctuations. They have exceptional dynamical proper¬ 
ties with exponential decay of spin correlation functions. 
Haldane predicted a unique many body singlet ground 
state and gapped magnetic excitations [3, 4]. The en¬ 
ergy gap between the singlet ground state and the 
excited state (triplet) is known as the Haldane gap. This 
is in contrast to a gapless continuum of multi-spinon ex¬ 
citations {S = 1/2) of its half-integer counterpart, the 
Spin-1/2 Heisenberg uniform AFM chain [3, 4]. The sin¬ 
glet ground state of the Haldane chain can be visualized 
in a simple way similar to the valence bond solid (VBS) 
state. Here, the VBS state is constructed by representing 
each S =1 spin as two separate S =1/2 spins which form 
a singlet with their counterparts on the two neighboring 
sites [5] as schematically shown in Fig. 1(a). Each ex¬ 
change bond carries exactly one valence bond, and the 
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periodicity of the underlying lattice remains unbroken. 

The excitation spectra of an isolated Haldane chain 
were investigated in great detail by several theoretical 
methods including quantum monte-carlo (QMC) [6], Ex¬ 
act diagonalization (ED) [7] and density matrix renor¬ 
malization group (DMRG) [8] as well as experimental 
techniques [9, 10]. The excitation spectra were well un¬ 
derstood as gapped dispersive excitation of magnons (S' 
= 1). The minimum gap of the excitation spectra (Hal¬ 
dane gap) appears at the AFM zone centre {qchain = tt)- 
DMRG calculations estimate the gap value Aq = 0.41J, 
where J is the intrachain antiferromagnetic coupling [8] . 
The one-magnon excitations disperse up to a maximum 
energy of ^ 2.4J at the zone boundary {qchain = 7J‘/2 
and 37r/2). A schematic excitation spectrum of an iso¬ 
lated isotropic Haldance chain is shown in Fig. 1(b). One 
of the interesting characteristics of the excitation spec¬ 
trum of Haldane chain is that the dispersion relation is 
asymmetric about qchain = and 37r/2 [Fig. 1(b) and 
(c)] since the translational symmetry of the lattice is not 
broken [9, 11]. This is in contrast to the conventional 
Neel antiferromagnets where the unit cell is doubled and 
qchain = 7r/2 and 37r/2 are of high symmetry points of 
the spin-wave dispersion relation. The intensity of the 
excitations [Fig. 1(b) and (c)] is well accounted by the 
single mode approximation [9] which shows a maximum 
at the zone center and decreases continuously away from 
q = TT and disappears rapidly beyond the zone boundary 



2 


(a) — 



Qchain /^ 


FIG. 1. (Color online) (a) A schematic representation of the 
singlet ground state of a Haldane chain in terms of a va¬ 
lence bond solid state, (b) A schematic representation (color 
plot) of the one-magnon excitation spectrum of an isolated 
Haldane chain over two Brillouin zones, (c) The dispersion 
relation (solid line) and the intensity profile (dashed line at 
the bottom) of the excitation spectrum in (b). 

[9, 10]. In addition to the one-magnon excitations, multi- 
magnon excitations are also present. Even- and odd- 
magnon excitations are predicted around the wavevec- 
tors qchain = 2mr, and mr, respectively, where n is an 
integer [8]. Thus the dynamics of the isolated Haldane 
chain are now well understood theoretically and experi¬ 
mentally. Our present research activities aim to explore 
and understand the influence of interchain interactions 
and anisotropy on Haldane chains. 

The presence of interchain interactions and anisotropy 
leads to complex behaviors and a rich phase diagram as 
mapped theoretically by Sakai and Takahashi [12]. The 
Haldane gap can even be suppressed by critical values 
of interchain interactions and anisotropy which results 
in three dimensional (3D) long-range magnetic ordering 
[13]. Experimental challenges are to investigate the ef¬ 
fects of interchain interactions and anisotropy on a model 
compound, to verify the theoretical predictions as well as 
to explore the critical properties. In this context, among 
the experimentally reported Haldane chain compounds, 
HNi2V208 where A = Sr or Pb are of particular interest 
since they have substantial interchain interactions and 
single ion uniaxial anisotropy [14, 15]. The compound 
SrNi2V208 is of present interest. The spin Hamiltonian 
for SrNi 2 V 208 can be defined as 

h = jY, S^.S^+l + J' Y, S^■S^+2 + JlY 

i i i,j 

+ ^D(5f)2,(l) 

i 

where the first and second terms are the nearest neighbor 
(NN) and next-nearest neighbor (NNN) intrachain inter¬ 
actions with exchange constants J, and J', respectively. 
The third term is due to interchain interactions where J_l 


is the effective interchain exchange constant. The fourth 
term is due to single-ion-anisotropy parameter D. 

SrNi 2 V 208 was proposed to be situated close to the 
phase boundary between the ordered and spin-liquid 
states in the Sakai-Takahashi phase diagram [14]. Based 
on initial powder inelastic neutron scattering (INS) stud¬ 
ies, Zheludev et al. [14, 15] proposed an ordered AEM 
ground state for SrNi 2 V 208 unlike the iso-structural 
compound PbNi2V208- However, a later nuclear mag¬ 
netic resonance (NMR) study on powder samples of 
SrNi2V208 by Pahari et al. [16] found a nonmagnetic 
spin-liquid (singlet) ground state. Moreover, our earlier 
bulk properties and preliminary neutron scattering mea¬ 
surements on a single crystal confirmed that SrNi 2 V 208 
has a non-magnetic spin-singlet ground state and gapped 
magnetic excitations [17]. 

Our previous single crystal study also revealed that 
a magnetic field induced quantum phase transition oc¬ 
curs from the gapped spin-liquid state to a 3D AEM 
ordered state [17]. The AEM ordered state is in con¬ 
trast to the theoretical prediction of a field-induced Lut- 
tinger liquid state for an isolated Haldane chain with¬ 
out anisotropy [18]. Moreover, when comparing the gap 
values to the critical magnetic fields, the Boson model 
[19, 20] correctly estimates gap value, in contrast to the 
Eermion and perturbative models [19, 21]. In contrast, 
for most of the studied Haldane chain compounds, which 
have planar anisotropy, the Fermion model estimates the 
gap values more accurately [22, 23]. In agreement with 
our finding, the gap values of the isostructural compound 
PbNi 2 V 208 which also has uniaxial anisotropy were also 
better described by the Bosonic model [24]. It was pre¬ 
dicted that the interchain interactions and the uniaxial 
single-ion anisotropy in SrNi 2 V 208 play a crucial role in 
determining the magnetic properties as well as in induc¬ 
ing the 3D AEM ordered state above the critical mag¬ 
netic field {He). However, the nature of interchain inter¬ 
actions and anisotropy, as well as their role in the spin 
correlations are still unknown and are of present research 
interest. 

In this paper, we report the results of single crystal in¬ 
elastic neutron scattering measurements on SrNi 2 V 208 
and DFT calculations to investigate the role of complex 
interchain interactions and single-ion anisotropy on the 
Haldane chain. Neutron scattering measurements were 
performed over a wide reciprocal space region to probe 
several Brillouin zones as well as over a wide energy range 
to map out the dispersions fully. The strengths of intra¬ 
chain and interchain interactions are estimated from the 
experimental excitation spectra by the Random phase 
approximation (RPA) method, and the exchange interac¬ 
tion geometry is proposed. The magnetic interactions are 
further investigated and verified by DFT based electronic 
band structure calculations. These calculations also pro¬ 
vide insight into the complex superexchange paths which 
give rise to these interactions. Finally, the position of the 
SrNi2V208 in the theoretical D-Jj^ (Sakai-Takahashi) 
phase diagram has been proposed showing where it lies 
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within the spin-liquid phase. 

II. METHODS 

Single crystals of SrNi 2 V 208 were grown in the Crys¬ 
tal Laboratory at Helmholtz Zentrum Berlin fiir Mate- 
rialen and Energie (HZB), Berlin, Germany. The start¬ 
ing material for crystal growth was prepared from high 
purity powders of SrCOa (99.994%, Alfa Aesar, Pu- 
ratronic), NiO (99.998%, Alfa Aesar, Puratronic) and 
V 2 O 5 (99.99% Alfa Aesar, Puratronic) by solid state re¬ 
actions. Crystal growth was performed in a four mir¬ 
ror type optical image furnace (Crystal Systems Corp., 
Japan) by the Traveling-solvent-floating-zone (TSFZ) 
technique. After growth the as-grown crystals were 
checked by x-ray and neutron Laue diffraction. Ground 
single crystals were also checked with x-ray powder 
diffraction (Brucker D 8 ) for phase purity. Detail of the 
crystal growth and characterization will be reported in a 
separate paper [25]. 

A powder neutron diffraction pattern of SrNi 2 V 20 g 
was recorded at 2 K by using the high resolution powder 
diffractometer E9 (A = 1.7982 A), at HZB, Germany. 
The diffraction pattern was analyzed by the Rietveld 
method using the FULLPROF software [26]. 

Excitation spectra were measured by single crystal in¬ 
elastic neutron scattering using the triple axis spectrom¬ 
eters with incident beams of thermal neutrons (PUMA 
at Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Ger¬ 
many and IN 8 at Institut Laue-Langevin (ILL), Greno¬ 
ble, France) as well as cold neutrons (PANDA at MLZ, 
Garching, Germany). For PUMA and INS, the measure¬ 
ments were performed with fixed final neutron wave vec¬ 
tor kf = 2.662 A and a Pyrolytic Graphite filter on 
the scattered side was used to remove higher order neu¬ 
trons. For PANDA, the final neutron wave vector was 
fixed to kf = 1.57 A and a cooled beryllium filter was 
used. Magnetic excitations were measured at base tem¬ 
peratures [at 3.1 and 3.5 K (using OCRs) for PANDA 
and PUMA, respectively, and at 1.5 K (using an orange 
cryostat) for IN 8 ]. All measurements were performed on 
a large cylindrical single crystal of mass ~ 2.5 g (diame¬ 
ter: ^ 6 mm and length: ~ 30 mm) which was mounted 
on an aluminum sample holder. 

Density functional theory (DFT) calculations were per¬ 
formed using the full-potential code fplo9. 07-41 [27]. 
For the generalized gradient approximation (GGA) and 
GGA-I-C/ calculations, we used the parameterization from 
Ref. 28. All calculations were performed for the experi¬ 
mentally determined crystal structure at 2 K. Nonmag¬ 
netic GGA calculations were done on a 16x16x16 k- 
mesh (648 points in the irreducible wedge). Wannier 
functions centered on Ni sites were computed accord¬ 
ing to the procedure described in Ref. 29. For the lo¬ 
cal X and y axes we used Ni~01 and Ni-04 vectors, 
respectively. Spin-polarized GGA-I-I/ calculations were 
done using symmetry-reduced cells (sp. gr. PI) on a 



FIG. 2. (Color online) The observed (filled circles) and 
calculated (solid black line) neutron diffraction patterns for 
SrNi2V208 measured at 2 K. The difference between observed 
and calculated patterns is shown by the thin blue line at the 
bottom. The vertical bars give the allowed Bragg peak posi¬ 
tions. 

2x2x2 fc-mesh. We used the on-site Coulomb repul¬ 
sion of Ud = 4.5 eV and the Bund’s exchange of Jd = 1 eV 
within the fully localized limit. 


III. RESULTS AND DISCUSSION 
A. Crystal structure 

First we discuss briefly the crystal structure of 
SrNi 2 V 208 which is important to understand the ex¬ 
change interaction geometry. Powder neutron diffrac¬ 
tion at 2 K using the E9 diffractometer at HZB, Berlin, 
Germany, was used to determine the low temperature 
structural parameters. The parameters were used for the 
analysis of the low temperature magnetic excitation spec¬ 
tra. The Rietveld refined diffraction pattern is shown in 
Fig. 2. The analysis confirms that the compound crystal¬ 
lizes in the tetragonal space group lAicd as found earlier 
at room temperature [30] . The refined lattice parameters 
are a = b = 12.135(1) A and c = 8.3159(1) A. The po¬ 
sitions of Sr, Ni and O were refined during the analysis. 
The refined values of their fractional coordinates at 2 K 
are given in Table 1. Since neutrons are not sensitive to 
the V nuclear position, the respective atomic coordinates 
were adopted from Ref. 31 and were kept fixed during the 
refinement. 

The spin chains in SrNi2V208 are formed by edge- 
shared NiOe (Ni^+; 3d®, S' = 1) octahedra along the c 
axis [Fig. 3(a)]. The screw chains are connected by non¬ 
magnetic VO 4 (V®+; 3d°, S = 0) tetrahedra. There are 
four screw chains per unit cell which are centered around 
(1/4, 1/4); (1/4, 3/4); (3/4, 1/4); and (3/4, 3/4) in the ab 
plane [Fig. 3(b)]. Each of these screw chains is 4-fold and 
contains four Ni ions along the c axis within one unit cell. 
The two diagonal chains rotate clockwise while the other 
two chains rotate anti-clockwise when propagating along 
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TABLE 1. The crystal structural parameters of SrNi2V20g 
determined by neutron diffraction at 2 K. The space group is 
74icd (110), a = 12.135(1) A, c = 8.3159(1) A. Full occupan¬ 
cies for all the atomic sites were considered. 


Atom 

Wyckhoff 

position 

xja 

via. 

zjc 

Biso (A=) 

Sr 

80 

0 

0 

0.25 

0.87(1) 

ya 

16b 

0.0802 

0.2596 

0.1485 

0.22 

Ni 

16b 

0.3304(4) 

0.1678(3) 

0.2793(1) 

0.40(4) 

01 

16b 

0.1476(5) 

0.0019(10) 

0.5051(15) 

0.49(1) 

02 

16b 

0.1601(9) 

0.3335(9) 

0.0214(11) 

0.77(3) 

03 

16b 

0.3196(7) 

0.1623(3) 

0.0371(1) 

0.26(3) 

04 

16b 

0.3326(6) 

0.0021(9) 0.3071(12) 

0.41(0) 


^ the values of the atomic coordinates and Bi^o were adopted 
from Ref. 31 



FIG. 3. (Color online) (a) The crystal structure of 
SrNi2V208. The Ni, V and O -ions are shown by green 
(large), yellow (medium) and red (small) spheres, respectively. 
The dotted lines indicate the dimensions of the unit cell. For 
clarity only two chains (out of four chains per unit cell) are 
shown, (b) Projection of the crystal structure on to the ab 
plane. Only the atoms and bonds are shown without polyhe- 
dra for clarity. Ni-O and V-O bonds are shown by green and 
yellow lines, respectively. Red arrows represent the direction 
of the rotation of the screw chains when propagating along 
the c axis. Neighboring chains rotate in opposite directions, 
(c) Arrangement of the Ni, V, and O atoms within a given 
ab plane (at c = 0). The direct distances between Ni ions 
along the diagonal directions are shown. The shortest diago¬ 
nal corresponds to the exchange interaction J4 which occurs 
via Ni-0..0-Ni pathway as indicated. 


the c axis. Furthermore, a phase shift of c/2 is present 
between diagonal chains [between (1/4, 1/4) and (3/4, 
3/4), and between (1/4, 3/4) and (3/4, 1/4)]. This re¬ 
sults in a 90° rotation of all the bonds within an ab plane 
[as shown in Fig. 3(c)] with respect to the neighboring ab 
plane along the c axis. This peculiar crystal structure 
provides multiple interchain interactions between neigh¬ 
boring chains [30]. The interchain interactions are not 
confined within the ab planes but can also have an out- 
of-plane component. For instance, the interchain inter¬ 
actions are possible between two Ni^+ ions (from two 
neighboring chains) which have an offset relative to each 
other of c/4 and c/2 along the c axis. In addition, next- 
nearest neighbor intrachain interactions, coupling two Ni 
ions separated by c/2, are also present due to the screw 
nature of the chains. 


B. Inelastic neutron scattering 

The excitation spectra [5'(g, w)] along the chain direc¬ 
tion, measured at 3.5 K using the PUMA spectrometer, 
show gapped magnetic excitations with strongly disper¬ 
sive modes up to ^ 23 meV [Figs. 4(a), 4(b) and 4(c)]. 
Complicated variations in the intensity of the excitation 
modes are found. In addition to the intensity modula¬ 
tion of Haldane chains (as discussed in Fig. 1), additional 
complexity arises due to the screw nature of the spin- 
chains in SrNi2V208. As discussed in the previous sec¬ 
tion the screw chains have four fold periodicity along the 
chain axis and each of the magnetic ions is shifted from 
the center of the chain axis. This leads to step-4 period¬ 
icity in the transverse displacement of the magnetic sites 
from the chain axis. Exact analytical calculation shows 
that the screw chains result in six modes with different 
structure factors [14]. All the individual modes have a 
periodicity of 4 r.l.u. along the chain direction (along the 
Qi) since there are four Ni ions per unit cell along the 
chain axis. However, the modes are shifted in the recip¬ 
rocal space (discussed later in detail). In a measurement 
along the Qi [such as along (1,0,1), (3,0,/), and (2,2,/) in 
Figs. 4(a)-(c)], the individual modes appear with a shift 
of integer units with respected to each other [Figs. 4(g)- 
(i)]. This leads to energy minima of one of the modes 
at each of the reciprocal lattice points (any integer com¬ 
binations of h, k, and /) which are refered to as AFM 
zone centers [Figs. 4(g)-(i)]. As the structure factor of 
individual modes is decided by the momentum transfer 
perpendicular to the chain axis (values of h and k) (dis¬ 
cussed later in detail) a substantial variation of relative 
intensity of the modes appears in the measured spectra 
[Fig. 4(a)-(c)]. For example, in the case of (3,0,/) fi¬ 
nite intensities are present only for the black [AFM zone 
center at (3,0,-3)] and the green [AFM zone center at 
(3,0,-1)] modes [Fig. 4(h)]. Intensities for the other four 
modes (red, blue, cyan, and magenta) are zero over the 
given reciprocal space. 

Three different energy minima (sizes ~ 2.3 meV, ^ 
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FIG. 4. (Color online) The measured and calculated magnetic excitation spectra. Left panel [(a), (b), and (c)] shows the 
measured magnetic excitation spectra along (1,0,Z), (3,0,Z), and (2,2,1), respectively, at 3.5 K using the triple axis spectrometer 
PUMA. The excitation spectra were obtained by combining several constant energy scans with an interval of 2 meV. Black 
dots represent the points of measurements. The panel at the centre [(d) (e), and (f)] shows the simulated magnetic excitation 
spectra along (1,0,Z), (3,0,/), and (2,2,1), respectively, as per the model discussed in the text and using the best fit parameters 
(Table II). The calculated spectra were convoluted with a Gaussian function and corrected for the magnetic form factor. The 
intensities are in arbitrary units, r.l.u. stands for reciprocal lattice units. Some additional intensities (sharp and diffuse) in the 
measured spectra are due to the contributions from phonon scattering. The right panel [(g), (h), and (i)] shows the calculated 
dispersion relations (solid curves) and intensities (dashed curves at the bottom) of each of the six individual modes that appear 
due to the screw nature of the chains (see text for details). For simplicity dispersion curves are shown without considering 
the anisotropy induced splittings. Different colors (blue, red, green, magenta, black and cyan) are used to represent different 
modes. 


3.6 meV and ~ 6.1 meV) are found at different AFM zone 
centers due to the complex interchain couplings. Fur¬ 
thermore, high resolution measurements using the cold 
neutron PANDA spectrometer [Figs. 5(a) and 5(c)] con¬ 
firm that each of the modes is split into two branches by 
anisotropy, giving the energy minima: (i) 1.57 ± O.OI and 
2.58 ± 0.01 meV, (ii) 3.07 ± 0.01 and 3.76 ± 0.01 meV, 
and (hi) 5.78 ± 0.01 and 6.18 ± 0.01 meV, respectively. 
The measured energy minima patterns for the AFM zone 
centers in two representative reciprocal planes (hkO) and 


(hkl) are shown in Figs. 6(a) and 6(b), respectively. For 
simplicity, only the average values of the anisotropy split 
gaps are presented. Complex patterns arise due to sev¬ 
eral competing interchain interactions (discussed later) 
and the screw nature of the chains. It is evident that all 
the energy minima values along some diagonal directions 
are the same (intermediate value of 3.6 meV) for partic¬ 
ular combination of h and k values i.e., even-odd for the 
I = 0 and both even-even and odd-odd for the / = 1. It 
may be worth to mention that the gap minima [red points 
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FIG. 5. (Color online) The measured and simulated magnetic 
excitation spectra perpendicular to the chain axis, (a) and (c) 
Measured excitation patterns using the PANDA spectrometer 
at 3 K, (d) using the INS spectrometer at 1.5 K, respectively. 
Black dots in (a) and (d) represent the points of measure¬ 
ments. The solid blue curves in (a) and (d) are the fitted 
dispersion curves, (b) and (e) The simulated spectra, corre¬ 
spond to (a) and (d), as per the model discussed in the text 
using the best fit parameters (Table II). For the comparison, 
simulated dispersion curves with J 4 = 0 (see text) are also 
shown in (b) and (e) by the dashed lines. The thin black lines 
through the measured points in (c) are the Gaussian fits to the 
data which are convoluted with instrumental resolution. The 
big solid circular points are the fitted peak positions and the 
thick lines are the htted dispersion curves as per the model 
discussed in the text. 


in Figs. 6(a) and 6(b)] become magnetic Bragg peaks in 
a magnetic field {H > He). 

The presence of dispersions [Figs. 5 and 6] perpendic¬ 
ular to the chain direction confirms that there are finite 
interchain interactions. The maximum bandwidth of the 
perpendicular dispersions (^ 4.0 meV) is about 20 % of 
the dispersion along the chain-direction. For the disper¬ 
sion along (/i,-5,l) [Figs. 5(d)], two modes (anti-modes) 
are present with the same periodicity (2 r.l.u.) in which 
one mode disperses between ~ 2.5-3.9 meV (lower energy 
mode) and the other mode disperses between ^ 3.9-6.6 
meV (higher energy mode). It may be noted that the 
bandwidth of the higher energy mode (^ 2.7 meV) is 
larger than that of the lower energy mode (^ 1.4 meV) 
which is found to be due to the competing interchain 
interactions (discussed later in detail). In addition, dis¬ 
persions are also found along the diagonal directions be¬ 
tween two neighboring reciprocal points having interme¬ 
diate energy minima values (mean value of 3.6 meV) as 
shown in Fig. 5(a), Fig. 6(c), and Fig. 6(d). These dis¬ 
persions indicate the presence of interchain interactions 
along the diagonal directions as found from the detailed 



(b)6 

5 

_ 4 
3 
2 
1 
0 



• • • • • 

#•••••# 

• • • • • • • 

*•••••• 


(h-h+5,0) [r.l.u.] 


FIG. 6 . (Color online) The measured energy minima patterns 
in the (a) {hkO) and (b) (hkl) reciprocal planes, respectively. 
The big (red), medium (green) and small (blue) points cor¬ 
respond to the average values of lowest (2.3 meV), medium 
(3.6 meV) and highest (6.1 meV) energy minima, respectively. 
The arrows indicate the measured dispersion along the diago¬ 
nal directions, (c) and (d) The measured dispersions between 
(1,1,1) and (2,2,1), and between (2,3,0) and (3,2,0) reciprocal 
points, respectively. The solid lines are the fitted dispersion 
curves. 


analysis (discussed later). 

To explain the observed features in the excitation spec¬ 
tra and to investigate the interactions and anisotropy, a 
model of the magnetic dynamical structure factor S{q, oj) 
(Eq. 6) was derived by the RPA method for the Hamil¬ 
tonian (Eq. 1) according to the method described by 
Zheludev et. al. in Ref. 14 and Ref. 15. The S{q, oj) for 
SrNi 2 V 208 was calculated in two steps; (i) first the calcu¬ 
lation of S {q, uj) for an equivalent Bravais spin-lattice of 
straight chains with an appropriate interaction topology 
[Fig. 7] and then (ii) the adaptation of the results for the 
more complex real crystal structure of SrNi 2 V 208 having 
screw chains with a fourfold periodicity [14]. 

The dynamical structure factors of a straight chain (in 
the single mode approximation) along the longitudinal 
(^ii) and transverse (S'j_) directions are given as 

( 2 ) 

where ly (= 2.49 J) is the spin-wave velocity, and Z = 1.26 
[14] and q is the reduced momentum transfer along the 
chain direction. For SrNi2V208, q = Qd^ = where 
Q (= ha* -\- kb* -|- Ic*) is the actual wave vector transfer, 
since there are four equivalent Ni^’*' ions per unit cell 
along the chain direction (c axis). 
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Intrachain interactions: 

J'. i/Ni..Ni = 2.866 A 
J'-. rfNi..Ni = 5.008 A 


Interchain interactions: 

J\'. rfNi..Ni = 6316 A 
6.387 A 

J 2 '. </Ni..Ni = 4.998 A 
J}’. dni..N\ = 5.851 A 
Ja: dm..Ni = 5.790 A 


FIG. 7. (Color online) The magnetic interaction topology 
of SrNi 2 V 208 (mapped onto equivalent straight chains). J: 
nearest neighbor intrachain interaction, J : next nearest 
neighbor intrachain interaction, Ji: between two Ni ions from 
neighboring chains having same c values (in the ab plane), J 2 : 
between two Ni ions from neighboring chains having offset 
of c/4, J 3 : between Ni ions from neighboring chains having 
offset of c/2, and J 4 : between Ni ions from diagonal chains 
within the ab plane. The Ni..Ni direct distances in the real 
crystal structure for corresponding interactions are also given. 


The dispersion relations for a coupled Haldane chain, 
like SrNi 2 V 2 08, can be written (within the RPA) as 

[?iw|| _L(q)]^ = + v'^si'n?{111/2) + cos^{nl/4) 

+ -^ J cos{itI )[1 — cos{7rl/2)] H — ■^J{Q)[^ — cos{iil/2)\ 

( 3 ) 

where A||_^ are the longitudinal and transverse Haldane 
gaps, respectively, for noninteracting chains. Two en¬ 
ergy gaps appear due to the zero-field splitting of the 
triplet states by single-ion anisotropy [32] . The relations 
between the energy gaps and single-ion anisotropy were 
proposed as [21, 33] 

Ax = Ao - 0.57TI and A|| = Ao -b 1.41Z9, (4) 

where Aq is the intrinsic Haldane gap (ss 0.41J) [ 8 ] in 
the absence of single-ion anisotropy. The 2nd and 3rd 
terms give the intrachain dispersions. The a (=1.1J) in 
Eq. 3 characterizes the asymmetry of the dispersion rela¬ 
tion [11, 15]. The asymmetry arises due to the fact that 
the translational symmetry is not broken in the Haldane 
ground state [ 6 ] . The fourth term is due to the contribu¬ 
tion from the NNN intra-chain exchange interaction J 
which is possible due to the screw chain crystal structure 
of edge sharing NiOe tetrahedra [Fig. 3]. The J{Q) in 
Eq. 3 is the combined Fourier transform of all interchain 
interactions between the magnetic sites in a crystallo¬ 
graphic unit cell (Bravais spin-lattice) [earlier referred as 
Jx in Eq. 1]. 

Now we consider the interchain interactions for the 
straight-chain spin-lattice. Only the interactions which 


TABLE II. The fitted values of the exchange interactions from 
neutron scattering. All values are in meV units. 


Interac¬ 

tion 

J 

y 

(2Jl- J2-bJ3) 

J 4 

Value 

D 

8.70T0.05 

-0.32±0.01 

0.149T0.004 

0.286T0.006 

0.32T0.01 


are possible in the real screw-chain crystal structure of 
SrNi 2 V 208 are considered (Fig. 7). The interactions are 
(i) Ji: between two Ni ions from neighboring chains hav¬ 
ing the same c values (in the ab plane) along the crystal¬ 
lographic a/b directions [dNi..Ni = 6.378/6.387 A in the 
real crystal structure and the coordination number {z) 
is four], (ii) J 2 '. between two Ni ions from neighboring 
chains offset by c/4 [(iNi..Ni = 4.998 A in the real crystal 
structure and the coordination number z = 2], (hi) J 3 : 
between two Ni ions from neighboring chains having an 
offset of c /2 [dNi..Ni = 5.851 A in the real crystal struc¬ 
ture and the coordination number z = 2], and (iv) J 4 : 
between two Ni ions from diagonal chains within the ab 
plane [(iNi..Ni = 5-'79 A in the real crystal structure and 
the coordination number z = 1]. Combining all these 
contributions, J{Q) becomes 

J{Q) =(2Ji + J2Cos{itI/2) + J 3 Cos{iTl))[cos{irh) -b cos( 7 rfc)] 
-b {J/i/2)cos{i:h)cos{iik) (5) 


Finally, the dynamical structure factors S{q, uj) for the 
real crystal structure of SrNi 2 V 208 are given by [Ref. 14] 


01 ^ 2 / 2 I s' fh l, , COS^/2lSiv?lll2 

^l±{q:^) = cos IpiCOS ’tp2S^^^j_{h,k,l) -\ - - - 

x[>S'|| j_(/l-b 1, k,l -b 1) -b -b 1, k,l + 3)] 

Sm^^/lCOS^'02 , ... r,' r, 1 

H-^- X [S^'^ j_{h,k -b 1 ,^ -b 1 ) -b 5||^_|_(h,fc -b 1,1 

-b3)] -b -b 1, fc -b 1, Z -b 2), ( 6 ) 

where 


and 1/2 = —k 
a a 


( 7 ) 


are the 3D structure factors that arise due to the screw 
nature of the spin chains and d (~ 0 . 8 a) is the offset of 
each Ni^+ ion along the a or b axis with respect to the 
central axis of the corresponding screw chains [(1/4, 1/4) 
or (1/4, 3/4) or (3/4, 1/4) or (3/4, 3/4)]. The dynamical 
structure factors S{q, oj) were simulated by using Eq. 6 
and compared to the measured spectra. By adjusting 
the values of the interactions a good solution could be 
found which is given in Table H and the corresponding 
simulated excitation patterns are depicted in Fig. 4 and 
Fig. 5 alongside the measured patterns. 

The analysis of the neutron scattering data reveals that 
the strongest interaction is the nearest neighbor intra¬ 
chain AFM interaction J = 8.7 meV. All other interac¬ 
tions are quite weak compared to J, although, compara¬ 
ble in strength to each other. Note that in the proposed 
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D/J 

FIG. 8. (Color online) The D-J± phase diagram for the 
weakly coupled Haldane chains (Ref. 12). The position of 
SrNi 2 V 208 obtained from the present study is shown by red 
square and compared with previous report (Ref. 14). Posi¬ 
tions of the other experimentally studied Haldane chain com¬ 
pounds are also shown as per Ref. 14. 


(a) 


a.H 


‘ 0 


H ;b 



FIG. 9. (Color online) (a) A schematic representation of the 
diagonal exchange in SrNi 2 V 20 g by a combination of inplane 
diagonal interchain interaction J 4 (dotted black line) and the 
NNN intrachain interaction J (solid yellow line), (b) The 
NN and NNN intrachain interactions J and J , respectively, 
in a single screw chain. 


model it is not possible to determine the individual val¬ 
ues of Ji, J 2 and J 3 since these interactions result in 
virtually identical dispersion relations. Only the linear 
combination of these parameters, (2Ji — J 2 -I- J 3 ) (since 
the coordination numbers for Ji, J 2 and J 3 are 2 ; = 4, z 
= 2 , and z = 2 , respectively) can be determined and is 
found to be AFM in overall. The minus sign of J 2 ex¬ 
presses the fact that it competes with Ji and J 3 if it has 
the same sign (e.g., AFM) but reinforces if it has opposite 
sign (e.g., FM). The diagonal interchain interaction J 4 is 
AFM and frustrated. J 4 strongly influences the disper¬ 
sion relation perpendicular to the chain direction, it (i) 
makes the band widths asymmetric [Fig. 5(d)], (ii) intro¬ 
duces dispersion along the diagonal directions [Fig. 5(a) 
and Figs. 6 (c)-(d)], as well as (iii) reduces the overall 
gap values. The NNN intrachain exchange interaction 
J is also AFM and frustrated which increases the gap 
values. The resulting effect of all the interchain interac¬ 
tions leads to the three different energy minima ^ 2.3, 
^ 3.6, and ^ 6.1 meV [Figs. 4, 5, and 6 ] at the AFM 
zone centers as compared to the intrinsic Haldane gap 
of (Aq = 0.41J)=3.57 meV. The fitted value of the uni¬ 
axial single-ion anisotropy parameter D is found to be 
—0.32 ± 0.01 meV. The value of the anisotropy parame¬ 
ter D is in good agreement with the value D = -0.29 meV 
obtained from the electron spin resonance measurements 
[34]. 

In SrNi 2 V 2 08, the values of all the interactions are 
such that the mean value (3.56 meV) of the intermediate 
energy minima (3.07 and 3.76 meV) is equal to the theo¬ 
retically predicted intrinsic Haldane gap value of (0.41J 
= 0.41 X 8.7 meV) 3.57 meV [8]. The first availabil¬ 
ity of good quality large single crystal has enabled us 
to perform such detailed investigations which give accu¬ 


rate values of exchange interactions as well as single-ion 
anisotropy. The present investigation estimates the val¬ 
ues zJ^jJ = 0.028 and D/J = —0.37 meV, respectively, 
and allows us to correctly position SrNi 2 V 208 in the the¬ 
oretical D-J± phase diagram [12], showing where it lies 
within the spin-liquid phase [Fig. 8]. 

Now we focus on the diagonal interchain interaction J 4 
which is of special interest. The presence of a diagonal in¬ 
teraction is strongly evident from the experimentally ob¬ 
served dispersions along the diagonal directions within 
the (hk) planes between two reciprocal lattice points 
having intermediate values (3.07 meV and 3.76 meV). 
For example, dispersions between (i) (2,3,0) and (3,4,0) 
[Fig. 5(a)], (ii) (1,1,1) and (2,2,1) [Fig. 6 (c)], as well as 
(iii) (2,3,0) and (3,2,0) [Fig. 6 (d)]. The simulated curves 
without considering the diagonal interchain interaction 
(i.e., J 4 = 0) reveal no dispersion between these recip¬ 
rocal lattice points. Such a simulated curve between the 
(2,3,0) and (3,4,0) reciprocal lattice points is shown by 
the dotted lines in Fig. 5(b). The presence of frustrated 
AFM J 4 is further evident from the fact that the band¬ 
width of the higher energy mode is larger than that of the 
lower energy mode for (h,-5,l) [Fig. 5(d)]. Here again, 
the patterns simulated with J 4 = 0 do not match with 
the observed excitation patterns. In contrast to the data, 
the simulation with J 4 = 0 results in a larger bandwidth 
of the lower energy mode than that of the higher energy 
mode as shown by the dotted curves in Fig. 5(e). 

Although the presence of diagonal interactions is con¬ 
firmed from the experimental data, the nature of the ex¬ 
change pathways is unclear for the real crystal structure 
of SrNi 2 V 208 . In the real crystal, the direct distances 
between Ni-ions along the diagonal directions are large 
~ 5.8 and 11.4, 8.9 and 9.1 A [Fig. 3(c)] which are far 
beyond the limit of any possible direct exchange inter- 
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actions. Moreover, along the diagonal directions, the 
Ni ions are not connected directly by any obvious path¬ 
ways (constructed from Ni-0 and V-0 bonds) [Fig. 3(c)] 
that can provide superexchange interactions. A superex¬ 
change interaction (J 4 ) is only possible between the Ni 
ions having shortest distance (5.8 A) via nontrivial Ni- 
0..0-Ni pathways [Fig. 3(c)]. This is confirmed from the 
DFT calculations (discussed later). 

The continuous propagation of this single J 4 interac¬ 
tion along a diagonal is not possible in SrNi 2 V 208 . If 
two Ni ions are separated along the diagonal by 5.8 A and 
coupled by J 4 , the next Ni ion in the same diagonal direc¬ 
tion will be situated far away (11.4 A) [Fig. 3(c)]. Thus 
there is no possible continuous exchange path along the 
inplane diagonal direction (1,1,0). Nevertheless, an effec¬ 
tive diagonal interaction along the body diagonal ( 1 , 1 , 1 ) 
direction may be plausible by a combination of J 4 and 
the NNN intrachain interaction (J ). Such an interaction 
scheme is shown in Fig. 9. As mentioned in the crystal 
structure section all the bonds within a ab plane rotate 
by 90 degrees with respect to the neighboring ab planes 
along the c axis due to the 4i glide symmetry of the space 
group lAicd. Therefore, the shortest diagonal distances 
(5.8 A) which are parallel to each other are shifted hy c/2 
along the c axis. These two paths can be connected to 
each other by the NNN intrachain interaction (J ) which 
occurs between two Ni ions shifted by c/2 along the c 
axis within a given chain. Thus, a combination of J 4 and 
J results in an effective diagonal interaction along the 
( 1 , 1 , 1 ) direction. 


C. DFT calculations 

To get further insight in the microscopic exchange in¬ 
teraction model, we employed DFT calculations. As dis¬ 
cussed in the previous section, the INS study can esti¬ 
mate the sums of different interchain couplings, yet in 
some cases their individual values remain unknown. The 
purpose of our DFT calculations is twofold. In DFT cal¬ 
culations, we can directly access all individual exchange 
couplings, our primary goal is to verify the exchange in¬ 
teractions provided by INS, to estimate their individual 
values, as well as supplement our microscopic magnetic 
model with the missing terms. Second, we are interested 
in the mechanism of the magnetic exchange, i.e. which 
orbitals (including those on nonmagnetic sites) have a 
sizable contribution to the magnetic exchange. 

Before we proceed to the DFT results, several remarks 
on the electronic structure of Ni^+ systems should be 
made. In SrNi 2 V 208 , the magnetic Ni atoms have a 
distorted octahedral environment. Hence, the electrical 
crystal field splits the d orbitals into two manifolds: t 2 g 
comprising xy, xz, and yz orbitals, and eg, comprising 
and 3z^ — orbitals. The d® electronic config¬ 
uration of Ni^+ implies that the lower-lying t 2 g orbitals 
are fully filled, while the eg orbitals are half-filled (one 
electron per orbital). Since the number of electrons is 



FIG. 10. (Color online) The valence band of SrNi2V208 as 
calculated using GGA. Total and atom-resolved density of 
states (DOS) are shown. The Fermi level is at zero energy. 
Inset: the orbital-resolved DOS for the Ni 3d states. 



r L P H r X 


/c-vector 

FIG. 11. (Color online) The GGA band structure of the Cg 
states (circles) and the fit using Ni Sd,^ 2 _y 2 and Sd^^ 2 _j. 2 - 
centered Wannier functions (dashed lines). The Fermi level is 
at zero energy. 


even, we expect a gap in the electronic excitation spec¬ 
trum (band insulator). This intuitive picture is fully sup¬ 
ported by our non-magnetic GGA calculations that yield 
a valence band dominated by Ni states split into two 
manifolds and a gap of ^^0.2 eV (Fig. 10). Such a small 
band gap renders SrNi 2 V 208 as a semiconductor, which 
is at odds with the yellow color of the crystals, indicat¬ 
ing the importance of strong electronic correlations that 
are severely underestimated in the GGA. We should also 
note that the valence band is dominated by the Ni states, 
yet there is a substantial contribution of O and V states 
especially in the eg channel (Fig. 10). Both O and V 
states play a crucial role for the intra-chain J exchange, 
as will be demonstrated below. 

The primitive unit cell of SrNi2V208, whose volume 
is twice smaller than the conventional body-centered 
tetragonal cell, contains 8 magnetic Ni atoms, i.e. we 
expect 8x5 = 40 bands with dominant Ni d character. 
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TABLE III. Leading transfer tij (in meV) and exchange Jij (in meV) integrals in SrNi2V208. The paths are indicated in 
Fig. 7. Interatomic Ni..Ni distances given in A correspond to the 2 K structure (Table I). Transfer integrals are provided for 
both magnetically active orbitals, 3z^ — and — y^. Note that due to the lack of inversion symmetry t^ 2 _y 2 ^^ 2 _y 2 A 
t^ 2 _y 2 ^^^ 2 _^ 2 . Terms with the absolute value larger than 5 meV are shown. The second last two columns show the GGA+17 
value for the magnetic exchange and the respective experimental values from INS. 


path 

dNL.Ni 

t3^2^2t^3^2^2 


tx2 —y2 r2 


J (GGA-G7) 

J (INS) 

interactions 

A 

2.866 

-94 

-111 

84 

216 

8.43 

8.7 

J 

A' 

5.008 

31 

35 

-34 

10 

0.21 

0.15 

J 

Ai 

A 2 

6.376 

12 



-6 

-0.11 


Ji 

J2 

6.387 

4.998 

14 

-8 

20 


-6 

7 

— 0.06 r/^ J 

-0.28 

-J2 + J3)=0.29] 

As 

5.851 

5 

-21 

10 

5 

-0.08 


Js 

A 4 

5.790 

-13 



45 

0.16 

0.32 

J4 



FIG. 12. (Color online) The Wannier function centered at 
the Ni orbital. The dominant exchange coupling J is 

realized by superexchange along the two inequivalent Ni-O- 
V-O-Ni paths indicated by arrows. 


Since the magnetism pertains to low-energy excitations, 
we restrict our analysis to the Cg states (Fig. 11) that 
lie in the vicinity of the Fermi energy. These 16 bands 
are projected onto a real-space basis of Ni-based Wannier 
functions, while the resulting overlap integrals tij param¬ 
eterize an effective two-orbital and 3z^—r‘^) tight- 

binding model. Going back to the A:-space, the tight- 
binding Hamiltonian yields excellent agreement with the 
GGA bands (Fig. 11). 

For each pair of Ni atoms, electron transfer can in¬ 
volve (i) two x^ — orbitals, (ii) two orbitals, 

or two orbitals of different character. Since SrNi2V208 
lacks inversion symmetry, the x'^ — y'^ —> and 

3z^ — x^ — y'^ couplings are not equivalent, hence 
there are four tij terms for each pair of Ni atoms i and 
j. The resulting transfer integrals ty are provided in 
Table III. All the leading terms correspond to the NN in- 
trachain coupling (J), and the x^ — y^ contribution to this 
coupling is by far the largest. The AFM J (~ 8.7 meV 
[Table II]) is surprisingly strong for an almost 90° Ni-0- 


Ni superexchange pathways between two intrachain NN 
Ni ions. According to the Goodenough-Kanamori rules 
[35, 36], a ferromagnetic interaction is preferable for 90° 
superexchange pathways. To pinpoint the cause of this 
dominance, we plot one of the Ni dy. 2 _y 2 -hased Wannier 
function in Fig. 12. By examining the nonlocal contribu¬ 
tions, we find that the bridging O atom in the Ni-O-Ni 
intrachain path only moderately contributes to the elec¬ 
tron transfer, while the main contribution is facilitated 
by the V d-orbitals and follows the nontrivial Ni-O-V- 
0-Ni superexchange pathways. 

Long-range superexchange paths involving empty d or¬ 
bitals are quite common in systems. For instance, 
Cu-O-V-O-Cu hopping gives rise to a sizable exchange 
of 330 K in the spin-chain compound AgCuV 04 [37]. 
A similar hopping process in / 3 -CU 2 V 2 O 7 leads to the 
honeycomb-lattice magnetism [38] , which is not observed 
in isostructural a-Gu 2 P 207 (coupled dimers. Ref. 39) and 
a-Gu 2 As 207 (alternating chains, Ref. 40). Hence, we 
can expect that the replacement of the d-element V in 
SrNi2V208 with a p element, such as P and As, would 
also drastically affect the magnetic coupling regime. 

After establishing the mechanism of the NN intrachain 
couplings J, we turn to the much smaller interchain terms 
as well as the NNN intrachain term. These transfer inte¬ 
gral terms are weak and amount to several meV at most, 
and this extremely low magnetic energy scale severely im¬ 
pedes the accuracy of DFT results. As follows from Ta¬ 
ble HI, several interchain couplings in SrNi2V208 have 
similar strengths. To estimate the numerical values of 
these small terms, we employ spin-polarized GGA-I-I/ 
calculations for magnetic supercells and map the result¬ 
ing total energies onto a classical Heisenberg model. The 
main drawback of this method is the dependence of the 
resulting exchange integrals on the Ud value used in 
GGA-I-C/. To eliminate this problem, we investigated a 
range of Ud values and evaluated the optimal Ud value 
of 4.5 eV, which reproduces the experimental value for J 
(8.7 meV) within the error bars. 

The resulting interchain exchange couplings are given 
in Table HI. A reasonably good agreement between the 
exchange interactions from the GGA-I-C/ and INS is 
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found. The possible intra and interchain interactions (for 
any transfer integral ty > 5 meV) agree well with the 
INS results and, hence, further long-range couplings in 
SrNi 2 V 208 can be neglected. In agreement with the INS 
results, the DFT calculations reveal the presence of the 
AFM NNN intrachain interaction J . The DFT calcula¬ 
tions provide the individual values of Ji, J 2 , and J 3 which 
are indistinguishable in the INS. In agreement with the 
INS, the overall sign of (2Ji — J 2 -|- J 3 ) is found to be 
antiferromagnetic. The DFT calculations further reveal 
that all the Ji, J 2 , and J 3 are ferromagnetic and the J 2 
is the strongest interaction among them. The presence of 
the antiferromagnetic J 4 is also confirmed between two 
Ni ions having the shortest distance of 5.8 A along the 
diagonal direction within a given ab plane. 


The exchange mechanism for the diagonal interaction 
J 4 was also revealed by the Ni (ia, 2 _y 2 -based Wannier 
function. By examining the nonlocal contributions, we 
find that the mechanism for J 4 is superexchange along 
Ni-0..0-Ni paths. It is quite efficient, because (i) the 
magnetically active x'^ — orbitals of the respective 
Ni atoms lie in the same plane and (ii) the direct 0..0 
distance (2.84 A) is short enough to facilitate overlap¬ 
ping of O 2p states. This is very similar to the situation 
in the edge-sharing chain cuprate LiVCu 04 [41], where 
the same mechanism gives rise to next-nearest-neighbor 
exchange interactions. 


IV. SUMMARY AND CONCLUSIONS 

In summary, our detailed single crystal inelastic neu¬ 
tron scattering study reveals that the ground state of the 
Haldane chain compound SrNi 2 V 20 g is a non-magnetic 
singlet, however, the low energy excitation spectra are 
significantly modified by the complex interchain interac¬ 
tions and single-ion anisotropy. The effective energy gap 
(1.57 meV) is reduced from the intrinsic Haldane gap 
value of 3.57 meV, and the system is pushed close to the 
phase boundary to long-range order. Comprehensive ex¬ 
perimental data along with DFT calculations reveal the 
important interactions, and estimate their strengths as 
well as the size of the single-ion anisotropy that allows 
SrNi 2 V 2 08 to be accurately positioned in the theoret¬ 
ical D-J± phase diagram. The position of SrNi 2 V 208 
within the spin-liquid phase is shown and compared with 
previous results. The presence of multiple competing 
inter-chain interactions is confirmed by both INS mea¬ 
surements and DFT calculations. The interchain inter¬ 
actions are quite weak (in the order of one tenth of meV) 
as compared to the strongest nearest neighbor intrachain 
interaction J = 8.7 meV. The DFT analysis reveals the 
mechanism of the exchange interactions, for example, (i) 
the strongest intrachain exchange J occurs via nontriv¬ 
ial Ni-O-V-O-Ni paths involving the empty d orbitals 
of the V ion and (ii) the diagonal interaction J 4 occurs 
via Ni-0..0-Ni paths with a direct overlapping of O 2p 
states. 
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